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Positrons from dark matter annihilation in the galactic halo: theoretical uncertainties 
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Indirect detection signals from dark matter annihilation are studied in the positron channel. We 
discuss in detail the positron propagation inside the galactic medium: we present novel solutions 
of the diffusion and propagation equations and we focus on the determination of the astrophysical 
uncertainties which affect the positron dark matter signal. We find dark matter scenarios and prop- 
agation models that nicely fit existing data on the positron fraction. Finally, we present predictions 
both on the positron fraction and on the flux for already running or planned space experiments, 
concluding that they have the potential to discriminate a possible signal from the background and, 
in some cases, to distinguish among different astrophysical propagation models. 

PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,95.30.Cq 



I. INTRODUCTION 



The quest for the identification of dark matter (DM), 
together with the comprehension of the nature of dark 
energy, is one of the most challenging problems in the 
understanding of the physical world. It is therefore of 
utmost importance to address the problem of the detec- 
tion of the astronomical DM with different techniques 
and in different channels: in underground laboratories, 
in neutrino telescopes, in large-area surface detectors as 
well as in space. Many efforts in both direct and indirect 
DM detection have been done in the last decade, and 
major breakthroughs are expected in the following years 
from the underground facilities and antimatter searches 
in space. In the same period, the LHC will provide crucial 
information on possible extensions of the Standard Model 
of particle physics, where the most viable DM candidates 
are predicted. We therefore are faced with the quest of 
signal predictions as detailed as possible, accompanied 
by a realistic estimation of their uncertainties. 
This paper deals with the indirect detection of DM 
through positrons from the DM pair annihilation inside 
the galactic halo. Secondary positrons and electrons are 
produced in the Galaxy from the collisions of cosmic- 
ray proton and helium nuclei on the interstellar medium 
[1] and are an important tool for the comprehension of 
cosmic-ray propagation. Data on the cosmic positron flux 
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(often reported in terms of the positron fraction) have 
been collected by several experiments [2, 3, 4, 5, 6, 7]. In 
particular, the HEAT data [2] mildly indicate a possible 
excess of the positron fraction (sec Eq. 32) for energies 
above 10 GeV and with respect to the available calcula- 
tions for the secondary component [1]. Different astro- 
physical contributions to the positron fraction in the 10 
GeV region have been explored [2], but only more accu- 
rate and energy extended data could shed light on the 
effective presence of a bump in the positron fraction and 
on its physical interpretation. Alternatively, it has been 
conjectured that the possible excess of positrons found 
in the HEAT data could be due to the presence of DM 
annihilation in the galactic halo [:•, : ]. This interpreta- 
tion, though very exciting, is at some point limited by 
the uncertainties in the halo structure and in the cosmic 
ray propagation modeling. Recently, it has been shown 
that the boost factor due to substructures in the DM 
halo depends on the positron energy and on the statisti- 
cal properties of the DM distribution [10]. In addition, 
it has been pointed out that its numerical values is quite 
modest [11]. 

The present work is about the issue of the propagation 
of primary positrons. We inspect the full solution of the 
diffusion equation in a two — zone model already tested 
on several stable and radioactive species [12] and quantify 
the uncertainties due to propagation models, in connec- 
tion with the positron production modes. Our results 
will be applied to experiments such as PAMELA and 
AMS-02, which are expected to bring a breakthrough in 
the cosmic antimatter searches and in the understanding 
of the positron component. In Sect. II we present the 
solutions to the diffusion equation with both the Green 
function formalism and the Bessel method, with a source 
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term due to the pair annihilations of DM particles. We 
introduce the diffusive halo function, the integral on the 
diffusive zone encoding the information relevant to cos- 
mic ray propagation through its fundamental parameters. 
In Sect. Ill we evaluate the uncertainties due to propa- 
gation models on the diffusive halo function, discussing 
the physical properties of the propagation parameter con- 
figurations giving the extremes of the uncertainty bands. 
The positron fluxes and the relevant positron fraction arc 
presented in Sect. IV, where we compare our results to 
existing data and elaborate predictions for present run- 
ning or planned experiments in space. In Sect. V we 
draw our conclusions. 

II. THE DIFFUSION EQUATION AND ITS 
SOLUTIONS 

The propagation of positrons in the galactic medium 
is governed by the transport equation 

^-y.{Ki^,E) V^}-±{biE)yj}^qi^,E) , (1) 

where iIj{x,E) denotes the positron number density per 
unit energy and q{K,E) is the positron source term. 
The transport through the magnetic turbulences is de- 
scribed by the space independent diffusion coefhcient 
K{x,E) = Koe^ where e = E/Eq and Eq = 1 GeV. 
Positrons lose energy through synchrotron radiation and 
inverse Compton scattering on the cosmic microwave 
background radiation and on the galactic starlight at a 
rate b{E) = Eq c'^/te where te = lO^^ s. The diffusive 
halo inside which cosmic rays propagate before escaping 



into the intergalactic medium is pictured as a flat cylinder 
with radius i?gai = 20 kpc and extends along the vertical 
direction from z ~ —L up to z = +L. The gaseous disk 
lies in the middle at z = and contains the interstel- 
lar material on which most of the cosmic ray spallations 
take place. The half-thickness L is not constrained by 
the measurements of the boron to carbon ratio cosmic 
ray fluxes B/C. Its value could be anywhere in the inter- 
val between 1 and 15 kpc. As cosmic rays escape from 
that diffusive zone (DZ) and become scarce in the inter- 
galactic medium, the density tp is generally assumed to 
vanish at the radial boundaries r = i?gai and z = ±L. 
Assuming steady state, the master equation (1) simplifies 
into 

Koe'A4^ + — A + q = , (2) 

oe [te J 

and may be solved by translating [s] the energy e into 
the pseudo-time 

i{e)^TE |t,(6) = ^| . (3) 

In this formalism, the energy losses experienced by 
positrons are described as an evolution in the pseudo- 
time t. As a consequence, the propagation relation (2) 
simplifies into the heat equation 

^ - i^oA^ = g(x,t) , (4) 

where the space and energy positron density is now ip = 
V whereas the positron production rate has become 
q = e^^^ q. 



In the Green function formalism, Ge+ (x, i? <— X5, Es) stands for the probability for a positron injected at X5 with 
the energy Es to reach the location x with the degraded energy E < Es, and the positron density is given by the 
convolution 

pEs=+oo /• 

V'(x,S)= / dEs d^xs G,+ (x,S^xs,^s) g(xs,Ss) . (5) 

JEs=E Jbz 



In the pseudo-time approach, the positron propagator 
may be expressed as 

Ge+(x,i?^Xs,i?5) = T?^G(x,t^X5,?s) , (6) 



where G is the Green function associated to the heat 
equation (4). Without any boundary condition, this heat 
propagator would be given by the 3D expression 




where t = t~-ts is the typical time including the diffusion 
process during which the positron energy decreases from 
Es to E. The distance between the source X5 and the 
observer x is Ax whereas the typical diffusion length as- 
sociated to f is Ad = V^KoT. In order to implement the 
vertical boundary conditions '0(±L) = 0, two approaches 
have been so far available. 

(i) In the regime where the diffusion length Ad is small 
with respect to the DZ half-thickness L, the method of 
the so-called electrical images consists in implementing 
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[8] an infinite series over the multiple reflections of the 
source as given by the vertical boundaries at +L and —L. 

(ii) In the opposite regime, a large number of images 
needs to be considered and the convergence of the scries 
is a problem. Fortunately, the diffusion equation along 
the vertical axis boils down to the Schrodinger equation 
- written in imaginary time - that accounts for the be- 
haviour of a particle inside an infinitely deep ID potential 
well that extends from z = —L to z = +L. The solution 
may be expanded as a series over the eigenstates of the 
corresponding Hamiltonian [10]. 

None of those methods deal with the radial boundaries 
at r = i?gai . The diffusive halo is here a mere infinite slab 
and not a flat cylinder. The Bessel approach which we 
present next remedies that problem and is an improve- 
ment with respect to the former Green formalism. 

A. The Bessel solution 



the cosmic ray density ^{r, z, e) as the Bessel series 



(8) 



Because the a^'s are the zeros of the Bessel function 
Jo, the cosmic ray density t/j systematically vanishes at 
the radial boundaries r = i?gai. The Bessel transforms 
Pi{z,e) fulfill the diffusion equation 



Kd^.P, -K^P, + —d, {e^P,} + Q,{z, e) = 0. (9) 

-"-gal '''E 



As the DZ is axisymmctric and since we will consider 
spherically symmetric source terms only, we may expand 



J 



The Bessel transform Qi of the source distribution q is given by the usual expression 

r-Rgai 



1 



Rl.i Jli^d Jo 



Jo(air'/i?gai) q{r, z, e) r dr. 



(10) 



r 



Each Bessel transform Pi(z,e) has to vanish at the The Bessel transform Pi(z, e) is Fourier expanded as 
boundaries z = —L and z = +L and may take any value 
in between. It can be therefore expanded as a Fourier 
series involving the basis of functions 



(13) 



(Pn{z) = sin(rifco z') 



(11) 



and the same expression holds for Qi{z,e) for which we 
where fco = 7r/2L and z' z + L. In our case, the need to calculate explicitly the Fourier coefficient 
DM distribution is symmetric with respect to the galactic 
plane and we can restrict ourselves to the functions ipn{z) 
with odd n ~ 2m + 1 



1 

;,ri(e) = J fn{z) Qi{z,e) dz 



(14) 



(pn{z) = (-1)™ cos(nfcoz) 



(12) 



J 



The Fourier transform of equation (9) involves the energy functions Pi,n(e) and Qi,n(e) 

1 



KnH^P.^n - K-j-P,,n + —d,{e^P,,n) + g»,„(z,e) = 



TE 



(15) 



At this stage, as for the Green approach, we can substitute the pseudo-time t for the energy e. By defining the new 
functions _Pi.„ = e^Pi^n and Qi,™ = £^~^Qi,n, we are led to the heat equation 



dt 



2l2 



^^al 



(16) 
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The solution to this ODE is straightforward 

hn{i)^ f Q,,„fe) cxp{-Q,„ (t-ts)} d?5 ■ (17) 



The argument of the exponential involves the diffusion 
length Ad through the pseudo-time difference f = t — 
as 




The cosmic ray positron density is given by the double 
expansion 

oo oo 

'^{'r,z,e) = Jo{a,r/Rgi,i) ifin{z) Pi,n{e-) , (19) 

i—1 n—1 

where 

P^,ni^) = ^£ Q,,„(es)exp{-Q,„ {i- ts)}des. (20) 

We eventually get the positron flux $e+ = 
(if.+ ipir, z, e)/AT: where the positron velocity {3^+ 
depends on the energy e. 

B. The source term for primary positrons 

Let us now consider the source term q (x, E) of the 
master equation (1). We are here interested in primary 
positrons, namely the ones that are produced by the pair 
annihilations of DM particles. According to the various 
supersymmetric theories, the annihilation of a DM pair 
leads either to the direct creation of an electron-positron 
pair or to the production of many species subsequently 
decaying into photons, neutrinos, hadrons and positrons. 
We have considered four possible annihilation channels 
which appear in any model of weakly interacting massive 
particles (WIMP). The first one is the direct production 
of a e^e~ pair and is actually generic for theories with 
extra-dimensions like the UED models [13, 14, 15]. The 
energy of the positron line corresponds to the mass of the 
DM species. We have alternatively considered annihila- 
tions into W+W~, r+r~ and bb pairs. These unstable 
particles decay and produce showers which may contain 
positrons with a continuous energy spectrum. Whichever 
the annihilation channel, the source term can be generi- 
cally written as 

q{^,E)=r^{av) (^j'/(,) . (21) 



The coefficient 77 is a quantum term which depends on 
the particle being or not self-conjugate : for instance, 
for a fermion it equals 1/2 or 1/4 depending on whether 
the WIMP is a Majorana or a Dirac particle. In what 
follows, we have considered a Majorana type species and 
taken t] = 1/2. The annihilation cross section is averaged 
over the momenta of the incoming DM particles to yield 
(av), the value of which depends on the specific SUSY 
model and is constrained by cosmology. We have actu- 
ally taken here a benchmark value of 2.1 x 10~^^ cm^ 
s^^ which leads to a relic abundance of fl^h^ ^ 0.14 (in 
agreement with the WMAP observations) under the hy- 
pothesis of dominant s-wave annihilation and by means 
of the relation: 

nh' = 8.5 • 10-11 GeV-^ ^ 3 ■ 10-^Ws-i 

9*s{xf) Xf-i (av) {av) 

(22) 

where Xf — m^/Tf ~ (20 -f- 25) with Tf the freeze-out 
temperature and where gn{xf) and g^^sixf) denote the 
effective number of degrees of freedom of the energy and 
entropy density at freeze-out, respectively. 

The DM mass is unknown. In the case of neu- 
tralinos, theoretical arguments as well as the LEP and 
WMAP results constrain this mass to range from a few 
GeV [16, 17, 18, 19] up to a few TeV. Keeping in mind 
the positron HEAT excess, we have chosen a neutralino 
mass of 100 GeV. We have also analyzed the positron sig- 
nal yielded by a significantly heavier DM particle of 500 
GeV. Finally, the energy distribution of the positrons 
produced in a single WIMP annihilation is denoted by 
/(e) = dNc+/dEc+ and has been evaluated with the help 
of the Pythia Montc-Carlo [20]. 

The only astronomical ingredient in the source 
term (21) is the DM distribution p(x) inside the Milky 
Way halo. We have considered the generic profile 

fro ^7 f l + (re/r.)" V^-^^/" 

where 7-0 = 8.5 kpc is the galactocentric distance of the 
solar system. Notice that r denotes here the radius in 
spherical coordinates. The solar neighborhood DM den- 
sity has been set equal to Pq = 0.3 GeV cm^''. Three 
profiles have been discussed in this work : an isothermal 
cored distribution [21] for which is the radius of the 
central core, the Navarro, Frenk and White profile [22] 
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Halo model 


a 


P 


7 


Ts [kpc] 


Cored isothermal [21] 


2 


2 





5 


Navarro, Frenk & White [22] 


1 


3 


1 


20 


Moore [2:!] 


1.5 


3 


1.3 


30 



TABLE I: Dark matter distribution profiles in the Milky Way. 



(hereafter NFW) and Moore's model [23]. The NFW and 
Moore profiles have been numerically established thanks 
to N-body simulations. In the case of the Moore profile, 
the index 7 lies between 1 and 1.5 and we have cho- 
sen a value of 1.3 - see Tab. I. The possible presence 
of DM substructures inside those smooth distributions 
enhances the annihilation signals by the so-called boost 
factor whose value is still open to debate. 

The positron flux at the Earth may be expressed as 



+ 00 



de5/(e5)/(AD)}, (24) 



where the information pertinent to particle physics has 
been factored out in 



K ~ 1] {(JV) 



(25) 



The diffusive halo integral / depends on the input energy 
ts and on the observed energy e through the diffusion 
length Ad given by 



Ad = ^KoTE 



1 



(26) 



In the Green formalism, the halo function / may be ex- 
pressed as the convolution of the reduced propagator G 
- see Eq. (6) and (7) - with the DM density squared 
{pI Pq)^ over the diffusive zone 



/(Ar 



d^xs G(0,e <- xs,es) 



DZ 



P(X5) 

Pq 



(27) 



Alternatively, in the Bessel approach, the halo integral / is given by the radial and vertical expansions 

00 00 

^(Ad)=^^ Jo{air/Rgi,i) Lpn{z) cxp j-Ci^n (f- ?s)| Ri,n , 

where the coefficients Ri „ arc the Bessel and Fourier transforms of the DM density squared {p/ . 

I 



(28) 



We insist again on the fact that the true argument of the 
halo function, whatever the approach followed to derive 
it, is the positron diffusion length Ad. This integral en- 
codes the information relevant to cosmic ray propagation 
through the height L of the diffusive zone, the normaliza- 
tion of the diffusion coefficient and its spectral index 
5. It is also the only relevant quantity concerning the 
DM distribution. The analysis of the various astrophys- 
ical uncertainties that may affect the positron signal of 
annihilating WIMPs will therefore be achieved by study- 
ing the behavior of /. 

C. The Bessel method versus the Green approach 

The diffusive halo integral /(Ad) may be calculated by 
using either the Bessel expansion method or the Green 
function approach. Wc investigate here the relevance of 
each as a function of the diffusion length Ad- 

To commence, the DM distribution is taken in Fig. 1 to 
be constant throughout the diffusive zone with p = pq. 
Both methods - Green or Bessel - do not give the same 
result as soon as L is large enough. Neglecting the radial 



boundary condition - the cosmic ray density vanishes at 
T = /?gai ~ leads to overestimate the halo function / when 
the diffusion slab is thick. This can be easily understood : 
if the slab is thin enough, a positron created near the 
radial boundary has a large probability to hit the vertical 
borders of the diffusive zone at ±L and hereby to escape 
into the intergalactic medium, never reaching the Earth. 
If so, the positron horizon does not reach radially the 
outskirts of the diffusive zone and the Green approach 
(long dashed curves) provides a very good approximation 
to the correct value of / as given by the Bessel expansion 
(solid lines). Conversely, if the slab is thick, we detect at 
the Earth a non negligible fraction of positrons produced 
near the radial boundary and the Green approximation 
is no longer acceptable. This is particularly true for the 
red curves of the L = 20 kpc case where the Green result 
largely overestimates the exact value. This justifies the 
use of the Bessel expansion method which improves upon 
the previous treatments of positron propagation and is 
one of the novelties of this article. 

As already discussed, a change in the normalization 
Kq or in the index 8 of the cosmic ray diffusion coeffi- 
cient leads only to a variation of the diffusion length Ad 
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1.2 - with radial boundaries (Bessel) 

without radial boundaries (Green) 




Diffusion length Xo=(4Kot)i/2 [kpc] 



FIG. 1: Influence of the radial boundary condition for a slab 
half-thickness L of 3, 10 and 20 kpc (i?gai = 20 kpc). The 
thicker the slab, the larger the error when neglecting the radial 
boundary. On the contrary, for small values of L, positrons 
produced near the radial outskirts of the diffusive halo es- 
cape into the intergalactic medium and do not contribute to 
the signal at the Earth. Implementing correctly the radial 
boundary condition is not relevant in that regime. 



through which those parameters appear. Notice that the 
relation that links the diffusion length Ad to the diffusive 
zone integral / is not affected by those modifications. On 
the contrary, the half-thickness L of the diffusive slab has 
a direct influence on the overall shape of / as a function 
of Ad as is clear in Fig. 2. In the left panel an isothermal 
distribution has been assumed whereas the right panel 
features the case of a NFW profile. For small values of L 
- see the green curve for which L = I kpc - the positron 
horizon is fairly limited. Because the positrons detected 
at the Earth merely originate from a very near region, 
the DM profile which we probe is essentially uniform. As 
in Fig. 1, the DZ integral / is unity below Ad ~ ^ and 
collapses for larger values of the diffusion length. For a 
thicker slab, the cosmic ray positron flux at the Earth 
gets sensitive to the center of the galaxy. That is why 
the halo integral / exhibits a maximum for a diffusion 
length ^5 — 7 kpc. a value close to the galactocentric 
distance tq = 8.5 kpc of the solar system. In both pan- 
els, the larger L, the more visible the bump. Notice also 
that the steeper the DM profile, the higher the maxi- 
mum. The curves featured in Fig. 2 point towards the 
importance of calculating correctly the influence of the 



DM located at the galactic center, when L is large. We 
therefore need to assess the relative merits of the Bessel 
and Green approaches in doing so. 

To achieve that goal, we have selected Moore's model 
with a very steep and dense DM central distribution. In 
the left panel of Fig. 3, the dashed curves are obtained by 
the Green method. The convolution (27) is numerically 
calculated by summing over the grid of elementary cells 
into which the diffusive halo has been split. The resolu- 
tion of that grid matters. For very small cells, the correct 
behavior of / is recovered and the dashed black curve is 
superimposed on the solid red line of the exact result de- 
rived with the Bessel expansion technique. However, the 
price to pay is an unacceptable CPU time. We had actu- 
ally to break the inner 1 kpc into 8 x 10^ cells in order for 
the integral (27) to converge. When the resolution of the 
grid is relaxed by increasing the size of the Green cells, 
the bump is dramatically underestimated. This is espe- 
cially clear for the dashed blue curve (labeled poor) where 
the DZ grid contains only a few 10'* cells. On the other 
hand, notice that even in that case, the correct result 
is obtained for a diffusion length smaller than ~ 3 kpc. 
In the right panel, we have concentrated on the Bessel 
method and numerically calculated the expansion (28). 
We have performed the summation up to a Bessel order 
of iVsessei and a Fourier order of iVharmonic ■ The exact 
result ~ featured by the solid red curve - incorporates a 
large number of modes and is once again obtained at the 
price of a long CPU time. If the expansion (28) is trun- 
cated earlier ~ see the dashed blue line - we observe that 
the correct value of / is completely missed when the dif- 
fusion length is small. In that regime, the positrons that 
are detected at the Earth originate mostly from the solar 
neighborhood. A large number A^Bessei of radial modes is 
needed in order for the Bessel transforms i?i_„ to inter- 
fere destructively with each other so that the influence 
of the galactic center is erased. For larger values of Ad, 
the exponential terms in Eq. (28) force the series over 
Ri,n to converge rapidly. The exact value of the halo 
integral / can be recovered even with as few terms as 
10 Bessel modes and 20 Fourier harmonics. Notice how 
well the peak at Ad ~ 7 kpc is reproduced by all the 
curves, whatever A^bossoI and iVharmonic- This peak was 
obtained with difficulty in the Green function approach. 
We therefore strongly advise to use the Green method as 
long as Ad is smaller than ^ 3 kpc whereas the Bessel ex- 
pansion technique should be preferred above that value. 
This prescription allows a fast and accurate evaluation of 
the halo function /. We can safely embark on an exten- 
sive scan of the cosmic ray propagation parameters and 
assess the theoretical uncertainties that may affect the 
positron DM signal at the Earth. 
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Diffusion length X^={4Kjy/'^ [kpc] 



Diffusion length Xo = (4K(,t)i/2 [kpc] 



FIG. 2: The halo convolution / is plotted as a function of the diffusion length Ad for various values of the slab half-thickness 
L. The left panel features the case of an isothermal DM distribution whereas a NFW profile has been assumed in the right 
panel - see Tab. I. When L is large enough for the positron horizon to reach the galactic center and its denser DM distribution, 
a maximum appears in the curves for Ad ~ J"©- 



D. The central divergence 



These conditions imply that po = p(ro) whereas 



Numerically derived DM profiles - NFW and Moore - 
exhibit a divergence at the center of the Milky Way. The 
density increases like r^^ for small radii - see Eq. (23) - 
but cannot exceed the critical value for which the WIMP 
annihilation timescale is comparable to the age of the 
galactic bulge. The saturation of the density typically 
occurs within a sphere of ~ 10~^ pc, a much shorter 
distance than the space increment in the numerical inte- 
grals, i.e., the Green convolution (5) and relations (10) 
and (14) for the Bessel expansion technique. Fortunately, 
this numerical difficulty can be eluded by noticing that 
the Green propagator Gg+ (0, e <— r ^ 0, eg) which con- 
nects the inner Galaxy to the Earth does not vary much 
over the central DM distribution. This led us to re- 
place inside a sphere of radius ro the r^'^ cusp with the 
smoother profile 

1 /2 

p*{r) ^ pq {1 + ai sinc(7ra;) -I- 02 sinc(27ra:)} , (29) 

where x = r/rp is the reduced radius. The coefficients po, 
ai and 02 are obtained by requiring that both the smooth 
density p* and its first derivative dp* /dr are continuous 
at tq. The other crucial condition is the conservation 
of the total number of annihilations within tq as the di- 
verging cusp p oc r"'' is replaced by the distribution p*. 



ai 



02 



02 = 87 



27 , 
TT^ - 9 + 67 

9(3 - 27) 



(30) 
(31) 



In Fig. 4, the halo integral / is plotted as a function of 
the diffusion length Ad in the case of the Moore profile 
and assuming a slab half-thickness L = 10 kpc. Within 
the radius tq, the central DM divergence has been re- 
placed cither by a plateau with constant density po or 
by the renormalized profile (29). In the case of the 
plateau, the maximum which / reaches for a diffusion 
length Ad ~ 7 kpc is underestimated even if values as 
small as 100 pc are assumed for tq. The larger that ra- 
dius, the fewer the annihilations taking place within tq 
as compared to the Moore cusp and the worse the miscal- 
culation of the halo integral. Getting the correct result 
featured by the solid red line would require a plateau 
radius so small that the CPU time would explode. On 
the contrary, we observe that the halo integral / is sta- 
ble with respect to a change of ro when the renormalized 
density p* is used. 
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1 10 1 10 



Diffusion length X^={4K^fy/^ [kpc] Diffusion length Xo = (4K(,f)i/2 [kpc] 

FIG. 3: The halo integral / is plotted as a function of the diffusion length Ad in the case of a Moore profile with L = 10 kpc. 
In the left panel, the results obtained with the Green function method are featured by the long-dashed curves and may be 
compared to the exact solution and its solid red line. In the right panel, the numbers A'sessei and A'harmonic of the eigenfunctions 
considered in the Bessel expansion (28) have been varied. The various curves reproduce astonishingly well the bump but diverge 
at small Ad when too few Bessel and Fourier terms are considered. 



III. PROPAGATION UNCERTAINTIES ON THE 
HALO INTEGRAL 

Following the prescription which has been given in the 
previous section, we can calculate accurately and quickly 
the halo integral / using cither the Green propagator 
method or the Bessel expansion technique according to 
the typical diffusion length Ad. We are now equipped 
with a rapid enough method for scanning the 1,600 
different cosmic ray propagation models that have been 
found compatible [12] with the B/C measurements. Each 
model is characterized by the half-thickness L of the dif- 
fusion zone and by the normalization Kq and spectral 
index S of the space diffusion coefficient. A large varia- 
tion in these parameters is found in [12] and yet they all 
lead to the same B/C ratio. The height L of the diffu- 
sion slab lies in the range from 1 to 15 kpc. Values of 
the spectral index 5 extend from 0.46 to 0.85 whereas the 
ratio Kq/L varies from lO"'^ to 8 x lO"'^ kpc Myr~^. 

In this section, we analyze the sensitivity of the 
positron halo integral / with respect to galactic propaga- 
tion. We would like eventually to gauge the astrophysi- 
cal uncertainties which may affect the predictions on the 
positron DM signal. A similar investigation - with only 
the propagation configurations that survive the B/C test 
- has already been carried out for secondary [25] and 



Model 


S 


Ko [kpc^Myr] 


L [kpc] 


MIN 


0.85 


0.0016 


1 


MED 


0.70 


0.0112 


4 


MAX 


0.46 


0.0765 


15 



TABLE II: Typical combinations of cosmic ray propagation 
parameters that are compatible with the B/C analysis [12] 
and which have been found [21] to correspond respectively to 
minimal, medium and maximal primary antiproton fluxes. 

primary [24] antiprotons. In the later case, three specific 
sets of parameters have been derived corresponding to 
minimal, medium and maximal antiproton fluxes - see 
Tab. II. 

Do these configurations play the same role for 
positrons ? Can wc single out a few propagation models 
which could be used later on to derive the minimal or 
the maximal positron flux without performing an entire 
scan over the parameter space ? These questions have 
not been addressed in the pioneering investigation of [9] 
where the cosmic ray propagation parameters have in- 
deed been varied but independently of each other and 
without any connection to the B/C ratio. 

In Fig. 5, we have set the positron detection energy E 
at a fixed value of 10 GeV and varied the injection energy 
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Moore DM halo profile (r^ = 30 kpc) 




1 10 



Diffusion length Xo=(4Kot)i/2 [kpc] 

FIG. 4: Same plot as before where the central DM profile 
within a radius ro is either a plateau at constant density po or 
the smooth distribution p* of Eq. (29). In the former case, the 
bump which / exhibits is significantly underestimated even 
for values of ro as small as 100 pc - solid dark blue - and 
drops as larger values are considered - solid light blue. On 
the contrary, if the DM cusp is replaced by the smooth profile 
p* , the halo integral no longer depends on the renormalization 
radius tq and the solid red and long-dashed black curves are 
superimposed on each other. 

Es from 10 GeV up to 1 TcV. The three panels corre- 
spond to the DM halo profiles of Tab. I. For each value 
of the injection energy Es, we have performed a com- 
plete scan over the 1,600 different configurations men- 
tioned above and have found the maximal and minimal 
values of the halo integral / with the corresponding sets 
of propagation parameters. In each panel, the result- 
ing uncertainty band corresponds to the yellow region 
extending between the two solid red lines. The lighter 
yellow domain is demarcated by the long-dashed black 
curves labeled MIN and MAX and has a smaller spread. 
The MED configuration is featured by the long-dashed 
blue line. In Fig. 6, the Moore profile has been chosen 
with four different values of the detection energy E. The 
corresponding uncertainty bands are coded with different 
colors and encompass each other as E increases. 

As Eg gets close to E, we observe that each uncertainty 
domain shrinks. In that regime, the diffusion length Ad 
is very small and the positron horizon probes only the 
solar neighborhood where the DM density is given by 
Pq. Hence the flagellate structure of Fig. 6 and a halo 
integral / of order unity whatever the propagation model. 



As is clear in Fig. 2, a small half-thickness L of the dif- 
fusion slab combined with a large diffusion length Ad im- 
plies a small positron halo integral /. The lower bound- 
aries of the various uncertainty bands in Fig. 5 and 6 
correspond therefore to parameter sets with L = 1 kpc. 
Large values of Ad are obtained when both the nor- 
malization Kq and the spectral index S are large - see 
Eq. (26). However both conditions cannot be satisfied to- 
gether once the B/C constraints are applied. For a large 
normalization Kq, only small values of S are allowed and 
vice versa. For small values of the detection energy E, 
the spectral index S has little influence on Ad and the 
configuration which minimizes the halo integral I corre- 
sponds to the large normalization Kq = 5.95 x 10^^ kpc^ 
Myr^^ and the rather small S = 0.55. For large values 
of E, the spectral index S becomes more important than 
Kq in the control of Ad. That is why in Fig. 6, the lower 
bound of the red uncertainty domain corresponds now to 
the small normalization Kq = 1.65 x 10"'^ kpc^ Myr~^ 
and the large spectral index 6 = 0.85. Notice that this 
set of parameters is very close to the MIN configuration 
of Tab. II. For intermediate values of E, the situation 
becomes more complex. We find in particular that for 
-E = 30 GeV, the halo integral / is minimal for the former 
set of parameters as long as Eg < 200 GeV and for the 
later set as soon as Eg > 230 GeV. In between, a third 
propagation model comes into play with the intermediate 
values Kq = 2.55 x 10"^ kpc^ Myr'^ and S = 0.75. It is 
not possible therefore to single out one particular com- 
bination of A'o and S which would lead to the minimal 
value of the halo integral and of the positron DM signal. 
The MIN configuration which appeared in the antiproton 
analysis has no equivalent for positrons. 

The same conclusion holds, even more strongly, in the 
case of the upper boundaries of the uncertainty bands. 
Whatever the DM halo profile, the panels of Fig. 2 feature 
a peak in the halo function / for large values of L and 
for a specific diffusion length Ag'*'^ ~ 7 kpc. At fixed E 
and Es, we anticipate that the maximal value for / will 
be reached for i = 15 kpc and for a diffusion length Ad 
as close as possible to the peak value Ad^'^- Two regimes 
can be considered at this stage. 

(i) To commence, the diffusion length Ad is below the 
critical value Ad'^'^ whenever the difference w(e) — v{es) 
is small enough - see the definitions (3) and (26). This 
condition is met in general when E and Eg are close to 
each other or when E is large. The largest possible value 
of Ad maximizes / and once again, we find two prop- 
agation models. For small E, the large normalization 
Kq = 7.65 X 10^2 kpc2 Myr-i is preferred with S = 0.46. 
We recognize the MAX configuration of Tab. II and un- 
derstand why the long-dashed black curves labeled MAX 
in the panels of Fig. 5 are superimposed on the solid 



10 



T. Delahaye, R. Lineros, N. Fornengo, F. Donato & P. Salati (2007) 



T. Delahaye, R. Lineros, N. Fornengo, F. Donato & P. Salati (2007) 



isothermal sphere (r^ = 5 kpc) 
E=10 GeV 



MAX 




a 

o 
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E=10 GeV 




Injection energy Eg [GeV] 
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Moore DM halo profile (r^ = 30 kpc) 




E=10 GeV 



Injection energy Eg [GeV] 

FIG. 5: In each panel, the halo integral I is plotted as a function of the positron injection energy Es whereas the energy E at 
the Earth is fixed at 10 GeV. The galactic DM halo profiles of Tab. I are featured. The curves labeled as MED correspond to the 
choice of cosmic ray propagation parameters which best-fit the B/G ratio [12]. The MAX and MIN configurations correspond 
to the cases which were identified to produce the maximal and minimal DM antiproton fiuxes [21], while the entire colored 
band corresponds to the complete set of propagation models compatible with the B/G analysis [12]. 



red upper boundaries. For large E, the spectral index 
6 dominates the diffusion length Ad and takes over the 
normalization Kq of the diffusion coefficient. The best 
model which maximizes / becomes then 6 = 0.75 and 
Ko = 2.175 X 10-2 kpc2 Myr-i. 



(ii) When the difference v{e) — v{es) is large enough, 
the diffusion length Ad may reach the critical value A^'*'^ 
for at least one propagation model which therefore max- 
imizes the halo integral. As E and Es are varied, the 
peak value of / is always reached when a scan through 
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FIG. 6; Same plot as before where the Moore DM profile has 
been selected. Four values of the positron detection energy E 
have been assumed. The flagella structure of this figure results 
from the widening of the uncertainty band as the detection 
energy E is decreased. 

the space of parameters is performed. This peak value 
corresponds to the maximum of the halo integral, hence 
a horizontal upper boundary for each of the uncertainty 
bands of Fig. 5 and 6. The set that leads to Ad = AgJ'*'^ 
is different for each combination of E and Es and is not 
unique. In the case of the NFW DM profile of Fig. 5, 
the halo integral / is maximized by more than 30 models 
above Es = 120 GeV. 

The complexity of this analysis confirms that the prop- 
agation configurations selected by B/C do not play the 
same role for primary antiprotons and positrons. The 
two species experience the propagation phenomena, and 
in particular energy losses, with different intensities. As 
pointed out in Ref. [2G] , the average distance traveled by 
a positron is sensibly lower than the one experienced by 
an antiproton produced in the halo. 



IV. POSITRON FLUXES 

Now that we have discussed in detail the solution of 
the propagation equation, and have identified and quan- 
tified the astrophysical uncertainties on the halo integral 
/, we are ready to apply our analysis to the theoreti- 
cal predictions for the positron signal at the Earth po- 
sition. The positron flux is obtained through Eq. (24). 



As stated in Sec. II B, we will not adopt specific DM 
candidates, but will instead discuss the signals arising 
from a DM particle which annihilates into a pure final 
state. We consider four different specific DM annihilation 
channels: direct e+e" production as well as W~^W~, bb 
and r+T^. The DM annihilation cross section is fixed 
at the value 2.1 x 10^^® cm'^ s^^ and we will con- 
sider the cases of a DM species with mass of 100 GeV 
and of 500 GeV. Generic DM candidates, for instance a 
neutralino or a sneutrino in supersymmetric models, or 
the lightest Kaluza-Klein particle in models with extra- 
dimensions, will entail annihilation processes with spe- 
cific branching ratios into one or more of these bench- 
mark cases. The positron flux in these more general sit- 
uations would simply be a superposition of the results 
for each specific annihilation channel, weighted by the 
relevant branching ratios and normalized by the actual 
annihilation cross section. 



Model 


S 


Ko [kpcVMyr] 


L [kpc] 


MED 


0.70 


0.0112 


4 


Ml 


0.46 


0.0765 


15 


M2 


0.55 


0.00595 


1 



TABLE III: Typical combinations of cosmic ray propagation 
parameters that are compatible with the B/C analysis [12]. 
The model MED has been borrowed from Tab. II. Models Ml 
and M2 respectively maximize and minimize the positron flux 
over some energy range - roughly above 10 GeV - the precise 
extent of which depends on the mass of the DM particle, on 
the annihilation channel and also on the DM profile. Note 
that Ml is the same as MAX in Tab. II but this is coincidental. 

In Fig. 7, the propagated positron flux <l>e+ - multi- 
plied by the square of the positron energy E for conve- 
nience - is featured as a function of E for a 100 GeV 
DM particle and a NFW density profile. The colored 
[yellow] area corresponds to the total uncertainty band 
arising from positron propagation. In all panels, it en- 
larges at low positron energy. This may be understood as 
a consequence of the behavior of the halo integral / which 
was analyzed in Sec. III. Positrons produced at energy 
Es and detected at energy E originate on average from 
a sphere whose radius is Ad. That positron sphere en- 
larges as E decreases and so does the uncertainty band. 
As positrons originate further from the Earth, the de- 
tails of galactic propagation become more important in 
the determination of the positron flux. On the contrary, 
high-energy positrons are produced locally and the halo 
integral / becomes unity whatever the astrophysical pa- 
rameters. Notice also that the uncertainty band can be 
sizeable and depends signiflcantly on the positron spec- 
trum at production. In the case of the e+e" line of 
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FIG. 7: Positron flux E'^^,+ versus the positron energy E, for a DM particle with a mass of 100 GeV and for a NFW profile 
- see Tab. 1. The four panels refer to different annihilation final states : direct e^e~ production (top left), bb (top right), 
W^W~ (bottom left) and r+r^ (bottom right). In each panel, the thick solid [red] curve refers to the best-fit choice (MED) 
of the astrophysical parameters. The upper [blue] and lower [green] thin solid lines correspond respectively to the astrophysical 
configurations which provide here the maximal (Ml) and minimal (M2) flux - though only for energies above a few GeV in the 
case of (Ml). The colored [yellow] area features the total uncertainty band arising from positron propagation. 



the upper left panel, the positron flux $e+ exhibits a 
strongly increasing uncertainty as E is decreased from 
down to 1 GeV. That uncertainty is one order of 
magnitude at E ~ 10 GeV, and becomes larger than 2 or- 
ders of magnitude below 1 GeV. Once again, the positron 
sphere argument may be invoked. At fixed detected en- 
ergy E, the radius Ad increases with the injected energy 
Es- We therefore anticipate a wider uncertainty band as 
the source spectrum gets harder. This trend is clearly 



present in the panels of Fig. 7. Actually direct produc- 
tion is affected by the largest uncertainty, followed by the 
T+T^ and W^W~ channels where a positron is produced 
either directly from the or from the leptonic decays. 
In the bb case, which is here representative of all quark 
channels, a softer spectrum is produced since positrons 
arise mostly from the decays of charged pions originat- 
ing from the quark hadronization. Most of the positrons 
have already a low energy Es at injection and since they 
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FIG. 8: Same plot as in Fig. 7 but with a DM particle mass of 500 GeV. 



are detected at an energy E Es, they tend to have 
been produced not too far from he Earth, hence a lesser 
dependence to the propagation uncertainties. The astro- 
physical configuration M2 - see Tab. Ill ~ provides the 
minimal positron flux. It corresponds to the lower bound- 
aries of the yellow uncertainty bands of Fig. 7. The Ml 
configuration maximizes the flux at high energies. For 
direct production and to a lesser extent for the t+t^ 
channel, that configuration docs not reproduce the up- 
per envelope of the uncertainty band in the low energy 
tail of the flux. As discussed in Sec. Ill, the response of 
<i>e+ to the propagation parameters depends on the de- 
tected energy E in such a way that the maximal value 



cannot be reached for a single astrophysical configura- 
tion. Finally, taking as a reference the median flux, the 
uncertainty bands extend more towards small values of 
the flux. In all channels, the maximal flux is typically a 
factor of ^ 1.5-2 times larger than the median predic- 
tion. The minimal flux features larger deviations with a 
factor of 5 for the bh channel aX E ^ \ GeV, of 10 for 
W+W- and of 30 for r+r". 

Fig. 8 is similar to Fig. 7 but with a heavier DM species 
of 500 GeV instead of 100 GeV. Since the mass is 
larger, so is on average the injected energy £'5. Notice 
that at fixed positron energy E at the Earth, the radius 
Ad of the positron sphere increases with Es- Wc there- 
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fore anticipate that the propagated fluxes are affected 
by larger uncertainties for heavy DM particles. Again, 
the maximal flux does not exceed twice the median flux, 
while the minimal configurations are significantly de- 
pressed. At the reference energy E = 1 GeV, reductions 
by a factor of 10 between the median and minimal pre- 
dictions are obtained for the bb channel and amount to 
a factor of 20 in the W^W~ case. They reach up to 2 
orders of magnitude for the direct positron production. 
In this large DM mass regime, the astrophysical configu- 
ration M2 does not reproduce by far the lower bound of 
the uncertainty band as it did for the 100 GeV case. The 
message is therefore twofold. 

(i) Once the positron spectrum at the source is chosen 
- and the corresponding branching ratios have been 
defined - the correct determination of the uncer- 
tainty which affects the flux at the Earth requires 
a full scan of the propagation parameter space for 
each energy E. The use of representative astro- 
physical configurations such as Ml and M2 would 
not provide the correct uncertainty over the entire 
range of positron energy E. 

(ii) However, specific predictions have to be performed 
for a given model of DM particle and a fixed set of 
astrophysical parameters. This is why fits to the ex- 
perimental data should be performed for each prop- 
agation configuration over the entire range of the 
measured positron energies E. The best fit should 
correspond to a unique set of astrophysical param- 
eters. This procedure is the only way to reproduce 
properly the correct and specific spectral shape of 
the flux. 

The effect induced by different DM profiles is pre- 
sented in Fig. 9, where the positron fluxes for the bb 
and W'^W~ channels are reproduced for the DM dis- 
tributions of Tab. I. The mass of the DM particle is 
fixed at = 100 GeV. Notice how steeper profiles en- 
tail larger uncertainties, especially for the upper bound. 
This is mostly due to the fact that for large values of 
L - for which larger fiuxes are obtained ~ the positron 
fiux is more sensitive to the central region of the Galaxy, 
where singular profiles like the NFW and Moore distri- 
butions have larger densities and therefore induce larger 
annihilation rates. On the contrary, the lower envelope 
of the uncertainty band is not affected by the variation 
of the halo profile. In this case, with typically small 
heights L, positrons reach the solar system from closer 
regions, where the three halo distributions arc very sim- 
ilar and do not allow to probe the central part of the 
Milky Way. 

Fig. 10 depicts the information on the positron flux un- 
certainty from a different perspective. The flux $e+ and 



its uncertainty band are now featured for fixed values of 
the detected energy E whereas the DM particle mass is 
now varied. The flux <I>g+ is actually rescaled by the prod- 
uct E^ rn^ <I>g+ for visual convenience. Each band corre- 
sponds to a specific detected energy E and consequently 
starts at = E. In the case of the W^W~ channel, 
the bands start at = mw because this channel is 
closed for DM masses below that threshold. The behav- 
ior of these bands can be understood from Fig. 6, where 
the halo function / is plotted for the same detected en- 
ergies, as a function of the injection energy Es- In the 
case of direct positron production, there is a simple link 
between the two figures, because the source spectrum in 
this case is just a line at Es = m~^. For the other 
channels the situation is more involved since we have a 
continuous injection spectrum with specific features as 
discussed above. The main information which can be 
withdrawn from Fig. 10 is that at fixed detection energy, 
the larger the DM mass, the larger the uncertainty. Let 
us take for instance a detection energy of i? = 3 GeV. For 
direct production, where Es = m^, increasing the DM 
mass translates into a larger radius Ad of the positron 
sphere. As a consequence, the uncertainty band enlarges 
for increasing masses. This occurs for all the annihilation 
channels, but is less pronounced for soft spectra as in the 
bb case. Similar conclusions hold for all the other values 
of E. 

Comparison with available data is presented in Fig. 11, 
12 and 13. In Fig. 11, the positron fraction 



TOT 



TOT 



(32) 



is plotted as a function of the positron energy E. The to- 
tal positron flux 'I'™"^ at the Earth encompasses the an- 
nihilation signal and a background component for which 
we use the results of Ref. [1] as parameterized in Ref. [S]) 

- see the thin solid [brown] lines. The electron flux is de- 
noted by <i>e- . The mass of the DM particle is 100 GeV 
and a NFW profile has been assumed. The data from 
HEAT [2], AMS [4, o], CAPRICE [C] and MASS [7] are 
indications of a possible excess of the positron fraction 
for energies above 10 GeV. Those measurements may be 
compared to the thick solid [red] line that corresponds 
to the MED configuration. In order to get a reason- 
able agreement between our results and the observations, 
the annihilation signal has been boosted by an energy- 
independent factor ranging from 10 to 50 as indicated in 
each panel. At the same time, the positron background 

- for which we do not have an error estimate yet - has 
been shifted upwards from its reference value of Ref. [8] 
by a small amount of 10%. As is clear in the upper left 
panel, the case of direct production offers a very good 
agreement with the potential HEAT excess. Notice how 



15 



T. Delahaye, R. Lineros, N. Fomengo, F. Donalo & P.Salad (2007) 




Positron energy [GeV] Positron energy [GeV] Positron energy [GeV] 



T. Delahaye, R. Lineios, N. Fomengo, F. Donato & P.Salad (20O7) 




Positron energy [GeV] Positron energy [GeV] Positron energy [GeV] 



FIG. 9: Positron flux + versus tlie positron energy E, for a DM particle mass of 100 GeV and for different halo density 
profiles : cored isothermal sphere [21] (left panels), NFW [22] (central panels) and Moore [2."i] (right panels) - see Tab. 1. The 
upper and lower rows correspond respectively to a bb and W^W' annihilation channel. In each panel, the thick solid [red] 
curve refers to the best-fit choice (MED) of the astrophysical parameters. The upper [blue] and lower [green] thin solid lines 
stand for the astrophysical configurations Ml and M2 of Tab. III. The colored [yellow] area indicates the total uncertainty 
band arising from positron propagation. 



well all the data points lie within the uncertainty band. 
A boost factor of 10 is enough to obtain an excellent 
agreement between the measurements and the median 
flux. A smaller value would be required for a flux at the 
upper envelope of the uncertainty band. The W~^W~ 
and T~^T~ channels may also reproduce reasonably well 
the observations, especially once the uncertainty is taken 
into account, but they need larger boost factors of the 
order of 30 to 40. On the contrary, softer production 
channels, like the bb case, are unable to match the fea- 
tures of the putative HEAT excess for this value of the 
DM particle mass. For all annihilation channels, the un- 
certainty bands get thinner at high energies for reasons 
explained above. They surprisingly tend to shrink also at 



low energies, a regime where the positron horizon is the 
furthest and where the details of galactic propagation are 
expected to be the most important. Actually, the annihi- 
lation signal turns out to be completely swamped in the 
positron background. In particular, the signal from direct 
production stands up over the background only for en- 
ergies larger than 5 GeV. The corresponding uncertainty 
on the positron fraction is at most of the order of 50% 
for energies between 10 and 20 GeV. In the other cases, 
the uncertainty bands are even thinner. Beware finally of 
the positron background which should also be affected by 
uncertainties due to secondary production processes and 
propagation. These uncertainties are not currently avail- 
able and there is clearly a need to estimate them in order 
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FIG. 10: For fixed values of the detected energy E, the uncertainty bands on the positron flux E^m'^^^+ are shown as a 
function of the mass of the DM particle. The energies considered in the figure are E = 3, 10, 30 and 100 GeV. Each band 
refers to one of those values and starts at — E. 



to properly shape theoretical predictions and to perform 
better study of the current and forthcoming data. Such 
an investigation would involve a comprehensive analysis 
and is out of the scope of the present article. 

Somehow different is the situation for larger masses 
of the DM candidate. Fig. 12 features the same infor- 
mation as Fig. 11, but now for = 500 GeV. In this 
case, all the annihilation channels manage to reproduce 
the experimental data, even the softest one bb. For direct 
production, the positron fraction is very large at energies 
above 40 GeV, where no data are currently available. 



This feature would be a very clear signature of DM anni- 
hilating directly into e^e~ pairs, with strong implications 
also on the nature of the DM candidate. For instance, 
bosonic dark matter would be strongly preferred, since 
Majorana fermionic DM, like the ncutralino, possesses a 
very depressed cross section into light fcrmions because 
of hclicity suppression in the non-rclativistic regime. As- 
trophysical uncertainties on the signal in this case show 
up more clearly than for the case of a lighter DM species, 
but still they are not very large. The drawback of having 
a heavier relic is that now the boost factors required to 
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FIG. 11: Positron fraction e'^ /{e~ +e^) versus the positron detection energy E. Notations are as in Fig. 7. In each panel, the 
thin [brown] solid line stands for the background [1, s] whereas the thick solid [red] curve refers to the total positron flux where 
the signal is calculated with the best-fit choice (MED) of the astrophysical parameters. Experimental data from HEAT [2], 
AMS [1, .-], CAPRICE [(.] and MASS [7] are also plotted. 



match the data are quite large. In Fig. 12 they range 
from 250 for the soft channel to 400 for the r+r~ case. 
Such large boost factors appear to be disfavored, on the 
basis of the recent analysis of Refs. [10, 11]. 

In Fig. 13, the positron flux (not the fraction) is com- 
pared to the available experimental data for a 500 GeV 
DM particle and a NFW profile. The solid thin [brown] 
line features the positron background which we shifted 
upwards by 10% with respect to the reference value of 
Rcf. [8]. The thick solid [red] line encompasses both 



that background and the annihilation signal which we 
calculated with the best-fit choice (MED) of the astro- 
physical parameters. Both curves have been derived as- 
suming solar modulation implemented through the force 
field approximation with a Fisk potential <f)p of 500 MV. 
The dashed [red] line instead corresponds to the total 
positron interstellar flux without solar modulation. No- 
tice that this curve is superimposed on the thick [red] 
line above ~ 10 GeV, a regime where cosmic ray propa- 
gation is no longer affected by the solar wind. A reason- 
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FIG. 12: Same plot as in Fig. 11 but with a mass of the DM particle of 500 GeV. 



ably good agreement between the theoretical predictions 
and the data is obtained, especially once the theoretical 
uncertainties on the annihilation signal are taken into ac- 
count. Notice that the spread of each uncertainty band is 
fairly limited as we already pointed out for the positron 
fraction. The reasons arc the same. 

Prospects for the future missions are shown in Fig. 14 
and 15. In Fig. 14, a 100 GeV DM particle and a NFW 
halo profile have been assumed. The median [red] curve 
corresponds to the prediction for the best-fit MED choice 
of astrophysical parameters whereas the upper [blue] and 
lower [green] lines correspond respectively to the Ml and 
M2 propagation models - sec Tab. III. Since wc arc deal- 



ing with predictions which will eventually be compared 
to the measurements performed over an entire range of 
positron energies, we have to choose specific sets of prop- 
agation parameters as discussed above in this Section. 
The upper and lower curves therefore do not represent 
the maximal uncertainty at each energy - though they 
may do so in some limited energy range - but instead 
they are "true" predictions for a specific set of propaga- 
tion parameters. Fig. 14 summarizes our estimate of the 
capabilities of the PAMELA detector [27] after 3 years 
of running. We only plotted statistical errors. We reach 
the remarkable conclusion that not only will PAMELA 
have the capability to disentangle the signal from the 
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background, but also to distinguish among different as- 
trophysical models, especially for hard spectra. Our con- 
clusion still holds for the bb soft spectrum for which the 
Ml, MED and M2 curves of the upper right panel differ 
one from each other by more than a few standard devi- 
ations. PAMELA could be able to select among them, 
even when systematical errors are included. 

In Fig. 15, the case of a 500 GcV DM particle is con- 
fronted with the sensitivity of AMS-02 for a 3-year flight. 
The possibility to disentangle the signal from the back- 
ground is also clearly manifest here, even once the as- 
trophysical uncertainties are included - provided though 



that boost factors of the order of 200 to 400 are possible. 
But, unless direct production is the dominant channel, a 
clear distinction among the various astrophysical models 
will be very difflcult because the Ml and M2 configura- 
tions are closer to the MED curve now than in the pre- 
vious case of a lighter DM species. Comparison between 
Figs. 14 and 15 clearly exhibits that at least below the 
TeV scale, the effect of the mass should not limit the 
capability of disentangling the annihilation signal from 
the background. More problematic is our potential to 
distinguish among different astrophysical models when 
the DM mass sizeably exceeds the 100 GeV scale. 
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FIG. 14: Predictions for PAMELA for a 3-year mission. Tlie positron fraction e^/{e~ + e"*") and its statistical uncertainty 
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V. CONCLUSIONS 

We have analyzed the positron signal from DM annihi- 
lation in the galactic halo, focusing our attention to the 
determination of the astrophysical uncertainties on the 
positron flux due to the positron propagation inside the 
galactic medium. 

Propagation of galactic cosmic rays has been treated 
in a two-zone model [12] and we have solved the diffusion 
equation for primary positrons both in the Green func- 



tion formalism and with the Bessel expansion method. 
We find that the most efficient way of dealing with 
positron propagation is to adopt the Green function 
method for values of the diffusion length Ad = \f^K^ 
smaller than ~ 3 kpc, and to employ the Bcssel function 
technique whenever Ad becomes larger. In this way the 
radial boundaries of the diffusion region (which are ne- 
glected in the Green function approach) can be properly 
coped with by the Bcssel expansion method. 

The propagation uncertainties on the halo integral 
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have been calculated for the ~ 1,600 different cosmic 
ray propagation models that have been found compatible 
[12] with the B/C measurements. These uncertainties are 
strongly dependent on the source and detection energies, 
Es and E. As Es gets close to E, wc observe that each 
uncertainty domain shrinks. In that regime, the diffusion 
length Ad is very small and the positron horizon probes 
only the solar neighborhood. In the opposite case, the 
uncertainty can be as large as one order of magnitude 
or even more. As positrons originate further from the 



Earth, the details of galactic propagation become more 
important in the determination of the positron flux. On 
the contrary, high-energy positrons are produced locally 
and the halo integral / becomes unity whatever the as- 
trophysical parameters. 

Inspecting directly the positron fluxes, typically, for a 
100 GeV DM particle annihilating into a bb pair, uncer- 
tainties due to propagation on the positron flux are one 
order of magnitude at 1 GeV and a factor of two at 10 
GeV and above. We find an increasing uncertainty for 
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harder source spectra, heavier DM, steeper profiles. 

The comparison with current data shows that the pos- 
sible HEAT excess is reproduced for DM annihilating 
mostly into gauge bosons or directly into a positron- 
electron pair, and the agreement is not limited by the as- 
trophysical uncertainties. A boost factor of 10 is enough 
to obtain an excellent agreement between the measure- 
ments and the median Hux, for a 100 GeV DM particle. 
A smaller value would be required for a flux at the upper 
envelope of the uncertainty band. 

We have finally drawn prospects for two interesting 3~ 
year flight space missions, like PAMELA, already in oper- 
ation, and the future AMS-02. We reach the remarkable 
conclusion that not only will PAMELA have the capabil- 
ity to disentangle the signal from the background, but it 
will also distinguish among different astrophysical mod- 
els, especially for hard spectra. For AMS-02 the possibil- 
ity to disentangle the signal from the background is also 
clearly manifest. We also wish to remind that improved 
experimental results on cosmic ray nuclei, expecially on 
the B/C ratio, will be instrumental to improve the deter- 
mination of the parameters of the propagation models, 
and will therefore lead to sharper theoretical predictions. 
This in turn will lead to a more refined comparison with 
the experimental data on the positron flux. Moreover, a 



good determination of the unstable/stable nuclei abun- 
dances like the ^'^Be/^Be ratio could shed some light on 
the local environment, which is certainly mostly relevant 
to the positrons. 

In the present paper we have thus presented the meth- 
ods and the practical tools to evaluate the primary 
positron fluxes in detailed propagation models. We have 
provided careful estimations of the underlying uncertain- 
ties and shown the extraordinary potentials of already 
running, or near to come, space detectors. 
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